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Methode de traitement d'images. 



DOMAINE TECHNIQUE 

invention concerne une methode de traitement d'images dans des images 
contenant des structures filiformes, comprenant en parall^le une etape de filtrage 
desdites images et une etape de decision destin6e ei selectionner les pixels d'images 
appartenant a une structure filiforme. L'Inventlon concerne egalement un programme 
d'ordinateur destine a la mise en oeuvre de cette methode ainsi qu'un appareil 
d'imagerie medicale comprenant un dispositif d'acquisition d'images et un systeme de 
traitement de ces images selon ladlte methode. L'inventlon peut etre appliqu6e, par 
exemple, dans I'industrie de I'imagerie medicale. 

ETAT DE LA TECHNIQUE ANTERIEURE 

De nombreuses images contiennent des structures filiformes 
interessantes. C'est le cas des images m^dicales, ou ces structures filiformes 
int6ressantes peuvent etre des vaisseaux sanguins ou un catheter. Le diam^tre 
typique d'un catheter varie entre 1 et 3 pixels et le diametre typique d'un vaisseau 
sanguin entre 5 et 20 pixels. 

Un systeme de traitement d'images sur une sequence d'images acquise 
par fluoroscopie aux rayons X est decrit dans le brevet US5809105. Ce systeme 
comprend un dispositif d'acquisition d'images numeriques, un dispositif de filtrage 
utilisant un filtre gaussien destine a 6limlrier le bruit provoque par les fluctuations 
al6atoires du signal electrique lors de I'acquisition de la sequence d'images et un 
dispositif de detection destine a distinguer les pixels appartenant a une structure 
interessante des pixels appartenant au fond d'image. Ce dispositif de detection est 
basee sur une methode de detection des contours utilisant le gradient et ne permet 
pas une bonne selection des structures filiformes interessantes, qui sont 
generalement difficiles a distinguer dans une image de fluoroscopie, car la puissance 
des rayons X utilises est faible et done I'image tr^s bruitee. 

EXPOSE DE L'INVENTION 

Un premier but de l'inventlon est de proposer une methode permettant 
d'ameliorer la selection des structures filiformes interessantes dans des images. 

Selon l'inventlon, une methode de traitement d'images telle que d^finie dans 
le paragraphe d'ouverture comprend les caracteristiques suivantes. L'^tape de 
decision comprend, en parallfele, une sous-6tape d'estimation de la direction de 
chaque pixel d'image ainsi qu'une sous-6tape d'analyse de connectivity de pixels 



voisins basee sur leurs directions a la sortie de la sous-etape d'estlmation de la 
direction de chaque pixel dimage et une sous-etape de selection de groupes de 
pixels, en fonction du resultat de la sous-etape d'analyse de connectivite de pixels 
voisins basde sur leurs directions, a la sortie de la sous-etape de filtrage. 

Un avantage de cette methode est qu'elle permet d'eiiminer les fausses 
alarmes tout en conservant les structures filiformes int6ressantes. Cette methode est 
robuste et produit des resultats stables sur de nombreux types d'images. Un autre 
avantage de cette methode est qu'elle permet de s6lectionner des structures 
filiformes epaisses, ce qui repr^sente un interet en imagerie medicale ou certaines 
structures comme des vaisseaux sanguins sont des structures filiformes epaisses. 

Un deuxieme but de Tinvention est de proposer un programme d'ordinateur 
et un systeme de traitement d'images destines a la mise en oeuvre de cette methode. 

Un troisieme but de I'invention est de proposer un appareil d'imagerie 
medicale comprenant un dispositif d'acquisition d'images et un dispositif permettant 
de mettre en oeuvre ladite methode. 

^invention prend les aspects suivants en consideration. La selection de 
structures filiformes int^ressantes dans des images est souvent difficile. C'est le cas 
dans les images medicales, qui peuvent etre acquises, par exemple, par fluoroscopie 
aux rayons X ou par resonance magnetique. 

Le principe suivant permet d'ameliorer la selection des structures 
filiformes interessantes dans une image. Une image numerique, constitute par un 
tableau de valeurs correspondant aux niveaux de gris de chaque pixel, est filtree 
grace a un filtre dont le role est de selectionner une grande partie des structures 
filiformes de I'lmage. La nature du filtre n'est pas critique, des lors qu'il remplit ce 
role. Ce filtre peut egalement effectuer un seuillage de Timage filtree, afin d'effectuer 
une premiere selection de pixels, dans le but de rtduire le nombre de pixels a traiter 
dans les Stapes suivantes. Ce filtrage joue a la fois le role du dispositif de filtrage et 
du dispositif de detection decrit en reference a Tart anterieur. 

Parallelement a ce filtrage, une estimation de direction est realisee sur 
chaque pixel de I'image numerique. Par exemple, on peut mesurer la direction 
tangentielle d'un pixel. En effet, les pixels voisins appartenant a une meme structure 
filiforme auront des directions tangentielles proches. Cette estimation de la direction 
de chaque pixel permet d'obtenir une image de direction, qui contient pour chaque 
pixel un numero de secteur angulaire correspondant a sa direction. Une analyse de 
connectivite de pixels voisins basee sur leurs directions est realisee sur I'image de 
direction, grace aux numeros de secteur angulaire contenus dans cette image de 
direction. Le but de cette analyse de connectivity de pixels voisins est d'obtenir une 
image connectee ou les pixels voisins ayant des numeros de secteur angulaire proches 
sont connectts entre eux. Une selection de groupes de pixels est ensuite appliquee, a 



partir de cette image connect^e, sur llmage filtree, afin d'eliminer ou de conserver 
certains groupes de pixels en fonction d'un parametre reglable. Par rapport h une 
detection de contour bas6e sur une methode de gradient, la m^thode d6crlte dans 
■'invention permet une selection plus fine des groupes de pixels appartenant h une 
structure interessante et ainsi d'eliminer une grande partie des fausses alarmes. En 
effet, alors qu'une methode de gradient permet de selectlonner tous les points d'une 
image prdsentant un certain contraste avec leurs plus proches voisins, sans tenir 
compte des voisins plus eloignes, la methode decrite dans I'invention permet de tenir 
compte des voisins plus 6loignes, ce qui est important dans la detection des 
structures filiformes, qui s'6tendent generalement sur une longueur de plusieurs 
dizalnes de pixels. 

On obtient ainsi une image dans laquelle les structures filiformes 
interessantes sont s6lectionnees. II est important de noter que les structures 
prdsentant un interet dans une image ne seront pas toujours les memes suivant 
I'application. Par exemple, dans une image medicale repr^sentant une structure 
arborescente de vaisseaux sanguins, le m6decin pourra avoir int^rSt h selectlonner le 
catheter, dont le diametre est inf^rleur aux vaisseaux principaux. Le fait d'utiliser des 
parametres reglables, 6ventuellement par le medecin, permet de selectlonner des 
structures filiformes de differentes epaisseurs, suivant I'application consld6r6e. 

L'invention et des caracteristlques additionnelles qui peuvent etre 
utilisees avec avantage pour mettre en oeuvre l'invention, seront d^crltes ci-dessous 
plus en detail par reference h des figures. 

BREVE DESCRIPTION DES FIGURES 

La figure 1 est un diagramme conceptuel illustrant les caracteristlques de 
l'invention ; 

La figure 2 illustre un exemple de filtre destine a selectlonner des structures 
filiformes dans une image ; 

Les figures 3a h 3d illustrent une technique basee sur I'utilisation du Hessien 
destin6e h I'estimation de la direction de chaque pixel d'image; 

La figure 4 illustre un exemple de dispositif d'analyse de connectivity de 
pixels voisins basee sur leurs directions ; 

Les figures 5a et 5b sont des diagrammes bloc illustrant respectivement un 
systeme de traitement d'images et un appareil d'Imagerle medicale destines h la mise 
en oeuvre de l'invention. 



MODE DE REALISATION DE L'INVENTION 

La figure 1 est un diagramme bloc illustrant les caracteristiques de 
rinvention. Une methode de traitement dimages comprend une 4tape de filtrage 
realis6e par un filtre 1, une sous-etape d'estimation de la direction de chaque pixel 
d'image r^alis^e par un dispositif d'estimation de direction 2, une sous-etape 
d'analyse de connectivity de pixels voisins realisee par un dispositif d'analyse de 
connectivity 3 et une sous-etape de selection de groupes de pixels realisee par un 
dispositif de selection 4. 

Une image numerisee NUMI est appliquee au filtre 1. Une image filtree 
RFI est obtenue, dans laquelle une grande partle des structures filiformes de I'image 
numerisee NUMI est selectionnee. Paraliyiement, I'image numerisee NUMI est 
appliquee au dispositif d'estimation de direction 2. Une image de direction DIRI est 
obtenue, qui contient un ensemble de num^ros de secteur angulaire correspondant 
aux directions de tous les pixels de I'image numerisee NUMI. Le dispositif d'analyse de 
connectivity 3 a pour entree I'image de direction DIRI. Le fonctionnement de ce 
dispositif d'analyse de connectivity 3 sera prycise en dytail sur la figure 3. Ce 
dispositif d'analyse de connectivity 3 permet d'obtenir une image connectee CONI 
dans laquelle les pixels de I'image de direction DIRI sont connect6s ou non en 
fonction de leurs numyros de secteur angulaire. L'image filtree RFI est ensuite 
appliquee au dispositif de selection de groupes de pixels 4, qui permet d'effectuer une 
syiection des pixels de I'image filtrye RFI appartenant h une structure filiforme 
intyressante, en fonction de I'image connectye CONI. On obtient ainsi une image 
traitye PRCI. 

La figure 2 illustre un exemple de filtre 1 destiny a syiectionner les 
structures filiformes dans une image. Dans I'exemple illustry sur cette figure, le filtre 

I s'applique sur un pixel P de I'image NUMI appartenant a une structure filiforme STR. 

II faut noter que le filtre 1 s'applique de la meme fagon sur toiis les pixels de I'image 
NUMI. En outre, dans la description des figures 2 a 4, les structures filiformes que I'on 
dysire syiectionner sont sombres par rapport au fond d'image. Un premier voisinage 
Nl de hauteur H et de longueur L est considyre d'un coty du pixel P. Un second 
voisinage N2 de memes dimensions est considyry de I'autre coty du pixel P. La taille 
de ces voisinages Nl et N2 est dyterminye grSce aux parametres ryglables H et L. 
L'ycartement GAP entre le premier voisinage Nl et le second voisinage N2, c'est a dire 
le nombre de pixels entre ces 2 voisinages, est choisi en fonction de la taille des 
structures filiformes STR que I'on dysire syiectionner. Plus I'epaisseur des structures 
filiformes STR que I'on dysire syiectionner est importante, plus I'ecartement GAP entre 
les deux voisinages Nl et N2 sera choisi important. Par exemple, si I'on dysire 
syiectionner des structures filiformes STR d'ypaisseur moyenne un pixel, ce qui 
pourrait etre le cas d'un catheter dans une image mydicale, on choisira un ecartement 



GAP entre les deux voisinages Nl et N2 de I'ordre de deux ou trois pixels. Si I'on 
desire s6lectionner des structures filiformes STR d'6paisseur moyenne comprises entre 
10 et 20 pixels, ce qui pourrait etre le cas d'une structure arborescente de vaisseaux 
sanguins dans une image m6dicale, on choisira un 6cartement GAP entre les deux 
voisinages Nl et N2 de I'ordre de 20 pixels. Le principe de ce filtre 1 est le sulvant. 
On effectue la somme Su des niveaux de gris de tous les pixels situes au dessus et en 
dessous du pixel P sur la hauteur H dans la structure STR. Pour chaque llgne de pixels 
de longueur L appartenant au voislnage Nl, symbolis^e par des pointilles, on calcule 
la moyenne Mu des niveaux de grIs de tous les pixels du voislnage Nl presents sur 
cette ligne. Pour chaque ligne de pixels de longueur L appartenant au voislnage N2, 
symbolisee par des pointings, on calcule la moyenne M12 des niveaux de gris de tous 
les pixels du voislnage N2 presents sur cette ligne. On calcule le minimum minide Mu 
et de Mi2. On effectue la somme Sv de tous les minimums mini des lignes de pixels de 
longueur L presentes sur la hauteur H. On calcule enfin la difference C entre 2u et Zv, 
afin d'obtenir le contraste au pixel P. Cette operation est effectuee pour tous les 
pixels de I'image NUMI, et on obtient I'image RFI dont la valeur de chaque pixel^^est 
egal au contraste C. 

L'avantage de ce filtre 1 est qu'il permet de s^lectionner des structures 
filiformes d'epaisseurs choisies et ainsi d'ellminer des structures plus 6paisses. Par 
exemple, en imagerie medicale, ce filtre permet de selectionner une structure 
arborescente de vaisseaux sanguins et d'attenuer des structures plus 6paisses comme 
une poche de sang. 

En outre, ce filtre peut effectuer une premiere selection permissive des 
pixels de I'image numeris6e NUMI appartenant a une structure filiforme int^ressante, 
sulvant un critere utilisant un parametre r6glable X. Par exemple, les pixels 
selectionn^s sont ceux dont la valeur est superieur i X fois la variance du bruit dans 
I'image filtree par le dispositif decrit ci-dessus. Le parametre X peut etre choisi plus 
ou moins important selon le degre de selection desire. Le role de cette premiere 
selection est de limiter le nombre de pixels pris en compte dans les etrapes 
ulterieures, afin d'eviter les problemes de surcharge de memoire de pile. 

Les figures 3a h 3d illustrent un exemple de technique destinee a 
I'estimation de la direction de chaque pixel d'image. La m^thode consiste a 
diagonaliser une matrice de Hessien HES, decrite sur la figure 3d, dont les elements 
diagonaux sont 6gaux a une premiere valeur Lxx et une deuxieme valeur Lyy et les 
autres elements sont egaux a une valeur Lxy. Pour chaque pixel de I'image numerisee 
NUMI, on calcule une valeur Lxx, une valeur Lyy et une valeur Lxy. La figure 3a 
illustre le calcul de la valeur Lxx. Pour un pixel donne de I'image numerisee NUMI, on 
considere un voislnage N55 de dimensions 5 pixels sur 5 pixels, centre sur ce pixel, et 



on effectue une convolution de ce voisinage N55 avec un masque Mxx. On obtient 
ainsi la valeur Lxx. La Hgure 3b lllustre le calcul de la valeur Lyy. Pour le pixel donne 
de nmage numeris6e NUMI, on effectue la convolution du voisinage N55 avec un 
masque Myy. On obtient ainsi la valeur Lyy. La figure 3c illustre le calcul de la valeur 
Lxy. Pour le pixel donne de I'image num^risee NUMI, on effectue la convolution du 
voisinage N55 avec un masque Mxy. On obtient ainsi la valeur Lxy. On effectue 
ensuite une diagonalisation de la matrice de Hessien HES et on retient le vecteur 
propre U associe a la plus petite valeur propre obtenue. Ce vecteur est le vecteur de 
direction tangentielle du pixel consid6re. On calcule ensuite la direction de ce vecteur, 
de telle sorte que I'angle obtenu est un angle compris entre 0 et 180 degres. On 
divise ensuite cet angle par 10 et on arrondit le resultat obtenu a rentier superieur. 
On obtient ainsi un numero de secteur angulaire K pour le pixel considere. Dans cet 
exemple, il existe 18 num^ros de secteur angulaire differents. II est possible d'utiliser 
un nombre different de num^ros de secteur angulaire, par exemple en divisant I'angle 
du vecteur U par 5 au lieu de 10. 

La figure 4 illustre un exemple de dispositif d'analyse de connectivite 3 
base sur la direction. L'exemple considere s'applique sur 6 pixels de I'image de 
direction DIRI numerot^s de 1 a 6 (P1,P6), possedant 6 vecteurs de directions 
associes numdrotes de 1 ^ 6 (U1,U6) et 6 numeros de secteur angulaire associes 
numerotes de 1 a 6 (K1,K6). Le principe de I'analyse de connectivite est le suivant. 
Pour un pixel donne PN (N etant un nombre entier compris ici entre 1 et 6), on 
considfere un voisinage NN s'^tendant dans la direction du vecteur UN, sur une 
longueur I r6glable. Dans le cas present, deux voisinages sont represent^, symbolises 
par des ellipses en pointings. Le voisinage Nl est un voisinage du pixel PI s'etendant 
dans la direction du vecteur Ul. Le voisinage N5 est un voisinage du pixel P5 
s'etendant dans la direction du vecteur U5. Supposons que le pixel PI soit le premier 
pixel de I'image de direction DIRI traits par le dispositif d'analyse de connectivite 3. 
La valeur de son numero de secteur angulaire Kl dans I'image de direction DIRI est 
comprise, par exemple, entre 0 et 18. Le pixel PI est alors labellise avec une 
etiquette egale a 19. Tous les pixels du voisinage Nl ayant un vecteur de direction 
associe proche de Ul seront labellises avec cette meme etiquette 19, on dira qu'ils 
appartiennent au label 19. Par exemple, on peut considerer qu'un vecteur UN d'un 
pixel PN du voisinage Nl est proche du vecteur Ul si les angles de ces 2 vecteurs 
different de moins de 10 degres, c'est a dire si le numero de secteur angulaire KN et 
le numero de secteur angulaire Kl different au maximum d'une unite. Dans I'exemple 
considere, les pixels P2 et P3 appartiennent au voisinage Nl. Le numero de secteur 
angulaire K2 etant 6gal au numero de secteur angulaire Kl et le numero de secteur 
angulaire K3 etant superieur de une unit6 au numero de secteur angulaire Kl, les 
pixels P2 et P3 appartiendront egalement au label 19. La meme operation est realisee 



sur des voisinages des pixels P2, P3 et P4. II est alors facile de constater que les 
pixels P2, P3, P4, et P5 appartiendront tous au label 19. Le pixel P6 poss^de un 
num^ro de secteur angulaire K6 sup6rieur de 6 unit& au numero de secteur angulaire 
K5. Par consequent, il appartiendra a un label different de 19, par exemple 20. 
L'operation est ainsi effectu^e de fagon recursive sur tous les pixels de I'image de 
direction DIRI, jusqu'^ ce que tous les pixels de I'image de direction DIRI soient 
labellis^s. On notera qu'une fois qu'un pixel de I'image de direction DIRI a ete 
labellis^, son etiquette ne peut plus changer lorsque I'analyse de connectivite se 
poursuit. On obtient ainsi I'image CONI, constitute des difftrents labels auxquels 
appartiennent les pixels. 

La derni^re etape de la mtthode consiste a selectionner des groupes 
de pixels dans I'image filtree RFI, afin de conserver les structures filiformes 
interessantes de cette image. Le principe de la selection de groupes de pixels est le 
suivant. On effectue les sommes ^ des contrastes (dans I'image filtree RFI) de tous 
les pixels appartenant au meme label dans I'image connectee CONL On realise cette 
operation pour tous les labels de I'image CONI, et on considfere le maximum Y de . 
tous les ^1 calcules. On ne retient que les pixels appartenant a un label tel que \i est 
superieur au quotient de Y par un parametre M rtglable. On obtient ainsi I'image 
traitee PRCI, constitute des pixels de I'image filtrte RFI qui ont it€ retenus apres la 
selection de groupes de pixels. Le choix du parametre M permet d'obtenir difftrents , 
degres de selection. En effet, si M est faible (en pratique de I'ordre de quelques 
unites), peu de groupes de pixels de I'image filtrte RFI seront retenus dans I'image 
traitee PRCI. Si M est grand (par exemple de I'ordre de quelques centaines eh 
pratique), la quasi-totalite des groupes de pixels de I'image filtree RFI sera retenue 
dans I'image traitee PRCI. II faut noter que pour des structures filiformes 6paisses 
comme une structure arborescente de vaisseaux sanguins, la selection de groupes de 
pixels donne de bons rtsultats du fait que les pixels voisins appartenant d une telle 
structure ont generalement la meme direction. 

La figure 5a illustre un systfeme de traitement d'images destine a la 
mise en oeuvre de I'invention. le systeme de traitement d'images comprend un 
dispositif 5 de traitement d'images selon I'invention, une station de travail 6, un ecran 
7, une souris 8 et un clavier 9, 

L'image numerisee NUMI est traitee par le dispositif 5 de traitement 
d'images, selon la methode decrite dans I'invention. Ce dispositif 5 de traitement 
d'images peut etre Implemente sous la forme d'un microprocesseur programmable 
faisant partie de la station de travail 6 ou d'un circuit comprenant des elements 
comme des filtres, des memoires ou des operateurs logiques destines a realiser les 
differentes etapes dtcrites dans I'invention. Uecran 7, la souris 8 et le clavier 9 



permettent la visualisation des images traitees et permettent le reglage des differents 
parametres reglables decrits dans Tinvention. Le systeme de traitement d'images peut 
etre relie h un dispositif de stockage afin de stocker les diff^rentes images obtenues 
au cours du traitement dimages. 

La figure 5b tllustre un appareil d'imagerie medicale destin^ a la mise 
en oeuvre de Tinvention. L'appareil comprend une source de rayons X 10, une table 11 
destinee a recevoir un patient 12 a examiner, un systeme optique 13 et un systeme 
de traitement d'images 14. 

Le systeme optique 13 permet d'acquerir des images numerlsees a 
partir des rayons X Issus de la source de rayons X 10 ayant traverse le patient 12. Le 
systeme de traitement d'images 14 peut etre, par exemple, le systeme de traitement 
d'images decrit sur la figure 5a. 

La description ci-dessus en reference aux figures illustre I'invention 
plutot qu'elle ne la limite. II est evident qu'il existe de nombreuses alternatives qui 
rentrent dans la portee des revendicatlons ci-jointes. A cet egard, quelques remarques 
sont faites ci-dessous. 

II existe de nombreuses manieres pour s^lectionner des structures 
filiformes dans une image. La figure 2 illustre seulement une mise en oeuvre possible 
dans laquelle le filtre 1 utilise deux voisinages d'un pixel, avec des dimensions 
reglables. II est possible d'utiliser d'autres fiitres pour selectionner des structures 
filiformes dans une image, prenant la forme de masques par exemple. 

II existe de nombreuses manieres pour effectuer une estimation de 
direction dans une image donnee. La figure 3 illustre seulement une mise en ceuvre 
possible dans laquelle I'estimation de la direction de chaque pixel s'appuie sur la 
diagonalisation d'une matrice de Hessien. II est possible d'utiliser d'autres techniques 
pour effectuer Testimation de la direction de chaque pixel, par exemple en utilisant 
une voie discrete. 

II existe de nombreuses manieres pour effectuer une analyse de 
connectivite de pixels voisins basee sur leurs directions dans une image donnee. La 
figure 4 illustre seulement une mise en oeuvre possible dans laquelle I'analyse de 
connectivity de pixels voisins bas^e sur leurs directions utilise un voisinage d'un pixel 
s'etendant dans la direction du pixel, la taille de ce voisinage etant reglable grace a 
un parametre I. II est possible d'utiliser d'autres techniques pour effectuer I'analyse 
de connectivite basee sur la direction, avec des voisinages differents. 



1. Methode de traitement dimages dans des images contenant des 
structures filiformes, comprenant en parall^le les etapes suivantes. 

une ^tape de filtrage desdites jmages ; 

une etape de decision destin^e a s^lectionner les pixels d'image appartenant 
^ une structure filiforme int6ressante, ladlte m6thode etant caracterisee en ce que 
I'etape de decision comprend en parallele une sous-etape d'estimation de la direction 
de chaque pixel d'image ainsi qu'une sous-etape d'analyse de connectivity de pixels 
voisins basee sur leurs directions ^ la sortie de la sous-etape d'estimation de la 
direction de chaque pixel d'image et une sous-etape de selection de groupes de 
pixels, en fonction du resultat de la sous-etape d'analyse de connectivity de pixels 
voisins basee sur leurs directions, a la sortie de la sous-etape de filtrage. 

2. M6thode de traitement d'images selon la revendication 1, caracterisee 
en ce que ladite etape de filtrage desdites images comprend une sous-etape de 
selection de pixels, les pixels d'une image selectionnes etant ceux dont le contraste 
est superieur a X fois la variance du bruit dans Timage, X ^tant un param^tre reglable 
par Tutilisateur. 

3. Methode de traitement d'images selon I'une des revendlcations 1 ou .2, 
caractyris^e en ce que ladite etape de filtrage desdites images utilise deux voisinages 
(Nl) et (N2) d'un pixel donn^, I'ecart (GAP) entre ces voisinages etant reglable par 
Tutilisateur. '^■y 

4. Methode de traitement d'images selon la revendication 3, caracterisee 
en ce que la hauteur (H) et la longueur (L) desdits voisinages sont reglables par 
Tutilisateur. 

5. Methode de traitement d'images selon Tune des revendlcations 1 a 4, 
caracterisee en ce que ladite sous-etape d'analyse de connectivity de pixels voisins 
basee sur leurs directions utilise un voisinage d'un pixel donne, ce voisinage 
s'etendant dans la direction du pixel consid^re, determinee lors de ladite sous-etape 
d'estimation de la direction de chaque pixel d'image. 

6. Methode de traitement d'images selon la revendication 5, caracterisee 
en ce que la longueur dudit voisinage est reglable par I'utilisateur. 

7. Methode de traitement d'images selon I'une des revendlcations 1 a 6, 
caracterisee en ce que ladite sous-etape de selection de groupes de pixels utilise un 
parametre M reglable par I'utilisateur, ce parametre M permettant de calculer la 
somme minimale des contrastes des pixels d'un groupe donne requise pour que ce 
groupe soit selectionne. 
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8. Programme d'ordinateur executable par I'intermediaire d'un 
processeur, destine a la mise en oeuvre d'une methode de traitement d'Images selon 
Tune des revendications 1 a 7. 

9. Systeme de traitement d'Images comprenant un ordinateur destine a 
5 la mise en oeuvre d'un programme d'ordlnateur selon la revendication 8 ou un circuit 

destine a la mise en oeuvre de la methode de traitement d'Images selon I'une des 
revendications 1 S 7, un dispositif de projection des images traitees selon ladlte 
methode et eventuellement un dispositif de stockage desdites images. 

10. Appareil d'imagerie medicate comprenant un dispositif d'acquisition 
10 d'Images et un systeme de traitement d'Images selon la revendication 9. 



Proc^de de traitement d'images. 



5 DOMAINE TECHNIQUE 

L'Invention concerne un precede de traitement dimages dans des 
images contenant des structures filifomnes, comprenant en parallele une etape de 
filtrage desdites images et une etape de decision destinee a seiectlonner les pixels 
d'images appartenant a une structure filiforme. invention concerne 4galement un 
10 programme d'ordinateur destine a la mise en oeuvre de ce proced^ ainsi qu'un 

appareil d'imagerie medicate comprenant un dfspositif d'acquisition dimages et un 
systeme de traitement de ces images selon ledit proc^de. L'lnventfon peut etre 
appliquee, par exemple, dans Industrie de I'Imagerfe medicale. 

15 ETAT DE LA TECHNIQUE ANTERIEURE 

De nombreuses Images contiennent des structures filiformes 
interessantes. C'est le cas des Images medicales, ou ces structures filiformes 
interessantes peuvent etre des vaisseaux sanguins ou un catheter. Le diametre 
typique d'un catheter varle entre 1 et 3 pixels et le diametre typlque d'un vaisseau 

20 sanguin entre 5 et 20 pixels. 

Un systeme de traitement d'images sur une sequence d'images 
acquise par fluoroscopie aux rayons X est decrit dans le brevet US5809105. Ce 
systeme comprend un dispositif d'acquisition d'images num^riques, un disposltif de 
filtrage utillsant un filtre gaussien destine a eliminer le bruit provoque par les 

25 fluctuations aleatoires du signal electrlque lors de I'acqulsition de la sequence 

d'images et un dispositif de detection destine a distlnguer les pixels appartenant a 
une structure interessante des pixels appartenant au fond damage. Ce dispositif de 
detection est basee sur une methode de detection des contours utillsant le gradient 
et ne permet pas une bonne selection des structures filiformes interessantes, qui 

30 sont generalement difficiles a distlnguer dans une image de fluoroscopie, car la 
puissance des rayons X utilises est faible et done I'image tres bruitee. 

EXPOSE DE LINVENTION 

Un premier but de Tinventlon est de proposer un procede permettant 
35 d'ameliorer la selection des structures filiformes interessantes dans des images. 

Selon I'Invention, un procede de traitement d'images tel que defini 
dans le paragraphe d'ouverture comprend les caracteristlques suivantes. L'etape de 
decision comprend, en parallele, une sous-etape d'estlmation de la direction de 
chaque pixel d'image aInsI qu'une sous-etape d'analyse de connecttvite de pixels 



voisins basee sur leurs directions a la sortie de la sous-etape d'estimation de la 
direction de chaque pixel d'Image et une sous-etap^e de selection de groupes de 
pixels, en fonction du resultat de la sous-etape d'analyse de connectivite de pixels 
voisins basee sur leurs directions, a la sortie de I'etape de flltrage. 
5 Un avantage de ce procede est qu'll permet d'eliminer les fausses 

alarmes tout en conservant les structures filiformes interessantes. Ce proced^ est 
robuste et prodult des resultats stables sur de nombreux types d'images. Un autre 
avantage de ce procede est qu'il permet de selectionner des structures filiformes 
epaisses, ce qui represente un interet en imagerie medicale ou certalnes structures 
10 comme des vaisseaux sanguins sont des structures filiformes epaisses. 

Un deuxieme but de Tinvention est de proposer un programme 
d'ordinateur et un systeme de traitement d'images destines a la mise en oeuvre de 
ce procede. 

Un troisieme but de I'inventlon est de proposer un apparell dimagerle 

15 medicale comprenant un disposltif d'acquisition damages et un disposltif permettant 
de mettre en oeuvre ledit procede. 

L'invention prend les aspects suivants en consideration. La selection 
de structures filiformes interessantes dans des images est souvent difficile. Cest le 
cas dans les images medicales, qui peuvent etre acquises, par exemple, par 

20 fluoroscopie aux rayons X ou par resonance magnetique. 

Le principe suivant permet d'ameliorer la selection des structures 
filiformes interessantes dans une image. Une image num^rique, constituee par un 
tableau de valeurs correspondant aux niveaux de gris de chaque pixel, est filtree 
grace a un filtre dont le role est de selectionner une grande partie des structures 

25 filiformes de I'lmage. La nature du filtre n'est pas critique, des lors qu'il remplit ce 
role. Ce filtre peut egalement effectuer un seuillage de I'image filtree, afin 
d'effectuer une premiere selection de pixels, dans le but de reduire le nombre de 
pixels a traiter dans les etapes suivantes. Ce flltrage joue a la fois le role du 
disposltif de flltrage et du disposltif de. detection decrit en reference a I'art 

30 anterieur. 

Paralletement a ce flltrage, une estimation de direction est reallsee 
sur chaque pixel de Tlmage numerlque. Par exemple, on peut mesurer la direction 
tangentielle d'un pixel. En effet, les pixels voisins appartenant a une mSme structure 
filiforme auront des directions tangentieiles proches. Cette estimation de la direction 
35 de chaque pixel permet d'obtenir une image de direction, qui contient pour chaque 
pixel un numero de secteur angulaire correspondant a sa direction. Une analyse de 
connectivite de pixels voisins basee sur leurs directions est reallsee sur Tlmage de 
direction, grace aux numeros de secteur angulaire contenus dans cette image de 
direction. Le but de cette analyse de connectivite de pixels voisins est d'obtenir une 
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image connectee ou les pixels voisins ayant des numeros de secteur angulaire 
proches sont connectes entre eux. Une selection de groupes de pixels est ensuite 
appliquee, a partir de cette image connectee, sur I'image filtree, afin d'eliminer ou 
de conserver certains groupes de pixels en fonction d'un parametre reglable. Par 
5 rapport a une detection de contour basee sur une methode de gradient, le procede 
decrit dans {'invention permet une selection plus fine des groupes de pixels 
appartenant a une structure interessante et ainsi d'eliminer une grande partie des 
fausses alarmes. En effet, alors qu'une methode de gradient permet de selectlonner 
tous les points d'une image pr^sentant un certain contraste avec leurs plus proches 

10 voisins, sans tenir compte des voisfns plus eloignes, le procede d^crlt dans 

invention permet de tenir compte des voisins plus eloignes, ce qui est important 
dans la detection des structures filiformes, qui s'etendent generalement sur une 
longueur de piusieurs dizalnes de pixels. 

On obtient ainsi une Image dans laquelle les structures filiformes 

15 interessantes sont selectionnees. II est important de noter que les structures 
presentant un Interet dans une image ne seront pas toujours les memes suivant 
Tapplication. Par exemple, dans une image medicale representant une structure, 
arborescente de valsseaux sanguins, le medecin pourra avoir interet a selectlonner 
le catheter, dont le diametre est inferieur aux valsseaux princlpaux. Le fait d'utilfser 

20 des parametres reglables, eventuellement par le medecin, permet de selectioriner 
des structures filiformes de differentes epaisseurs, suivant {'application consideree. 

LMnvention et des caracteristiques additionnelles qui peuvent etre 
utilisees avec avantage pour mettre en ceuvre Tinvention, seront decrites ci-dessous 
plus en detail par reference a des figures. 

25 

BREVE DESCRIPTION DES FIGURES 

La figure 1 est un diagramme conceptuel illustrant les caracteristiques 
de rinvention ; 

30 La figure 2 illustre un exemple de filtre destine a selectionner des 

structures filiformes dans une Image ; 

Les figures 3a a 3d iliustrent une technique basee sur I'utllisatlon du 
Hessien desttnee a Testlmation de la direction de chaque pixel d'image; 

La figure 4 illustre un exemple de dispositif d'analyse de connectivite 
35 de pixels voisins basee sur leurs directions ; 

Les figures 5a et 5b sont des diagrammes bloc illustrant 
respectivement un systeme de traitement d'Images et un appareil d'imagerie 
medicale destines a la mise en oeuvre de {'invention. 



MODE DE REAUSATION DE LINVENTION 

La figure 1 est un diagramme bloc illustrant les caracteristiques de 
5 rinvention. Un procede de traitement d'images comprend une etape de.filtrage 
realisee par un filtre 1, une sous-etape d'estlmation de la direction de chaque pixel 
damage realisee par un disposittf d'estlmation de direction 2, une sous-etape 
d'analyse de connectlvite de pixels voisins realisee par un dispositif d'anaiyse de 
connectivite 3 et une sous-etape de selection de groupes de pixels realisee par un 

10 dispositif de selection 4. 

Une image numerisee NUMI est appliquee au filtre 1. Une image 
filtree RFI est obtenue, dans laquelle une grande partle des structures filiformes de 
{'image numerisee NUMI est selectionnee. Parallelement, I'image numerisee NUMI 
est appliquee au dispositif d'estimation de direction 2. Une image de direction DIRI 

15 est obtenue, qui contient un ensemble de numeros de secteur angulaire 

correspondant aux directions de tous les pixels de I'image numerisee NUMI. Le 
dispositif d'anaiyse de connectivite 3 a pour entree I'image de direction DIRI. Le 
fonctionnement de ce dispositif d'anaiyse de connectivite 3 sera precise en detail sur 
la figure 3. Ce dispositif d'anaiyse de connectivity 3 permet d'obtenir une image 

20 connectee CONI dans laquelle les pixels de I'image de direction DIRI sont connectes 
ou non en fonction de leurs numeros de secteur angulaire. L'image filtree RFI est 
ensuite appliquee au dispositif de selection de groupes de pixels 4, qui permet 
d'effectuer une selection des pixels de I'image filtree RFI appartenant a une 
structure filiforme Interessante, en fonction de I'image connectee CONI. On obtient 

25 ainsi une Image traltee PRCI. 

La figure 2 illustre un exemple de filtre 1 destine a selectionner les 
structures filiformes dans une image. Dans I'exemple illustre sur cette figure, le filtre 
1 s'applique sur un pixel P de I'image NUMI appartenant a une structure filiforme 
STR. II faut noter que le filtre 1 s'applique de la meme fagon sur tous les pixels de 

30 I'image NUMI. En outre, dans la description des figures 2 a 4, les structures 

filiformes que I'on desire selectionner sont sombres par rapport au fond d'image. Un 
premier voisinage Nl de hauteur H et de longueur L est considere d'un cote du pixel 
P. Un second voisinage N2 de memes dimensions est considere de I'autre cote du 
pixel P. La taille de ces voisinages Nl et N2 est determinee grace aux parametres 

35 reglables H et L. L'ecartement GAP entre le premier voisinage Nl et le second 

voisinage N2, c'est a dire le nombre de pixels entre ces 2 voisinages, est choisi en 
fonction de la taille des structures filiformes STR que Ton desire selectionner. Plus 
I'epaisseur des structures filiformes STR que I'on desire selectionner est importante, 
plus l'ecartement GAP entre les deux voisinages Nl et N2 sera choisi important. Par 
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exemple, si I'on desire selectionner des structures filiformes STR d'epaisseur 
moyenne un pixel, ce qui pourralt etre le cas d'un catheter dans une Image 
medicale, on cholsira un ecartement GAP entre ies deux volsinages Nl et N2 de 
I'ordre de deux ou trois pixels. Si I'on desire selectionner des structures filiformes 
5 STR d'epaisseur moyenne comprises entre 10 et 20 pixels, ce qui pourralt etre le cas 
d'une structure arborescente de vaisseaux sanguins dans une image medicale, on 
choisira un ecartement GAP entre Ies deux volsinages Nl et N2 de Tordre de 20 
^ pixels. Le princIpe de ce flltre 1 est le suivant. On effectue la somme Zy des niveaux 
de grls de tous Ies pixels situes au dessus et en dessous du pixel P sur la hauteur H 

10 dans la structure STR. Pour chaque llgne de pixels de longueur L appartenant au 
volsmage Nl, symbolisee par des polntllles, on calcule la moyenne Mu des niveaux 
de grls de tous Ies pixels du voisinage Nl presents sur cette Hgne. Pour chaque ligne 
de pixels de longueur L appartenant au voisinage N2, symbolisee par des pointilles, 
on calcule la moyenne Mu des niveaux de gris de tous Ies pixels du voisinage N2 

15 presents sur cette llgne. On calcule le minimum minide Mu et de M^. On effectue la 
somme Sy de tous Ies minimums mini des lignes de pixels de longueur L pr^sentes 
sur la hauteur H. On calcule enfin la difference C entre et ILy, afin d'obtenir le 
contraste au pixel P. Cette operation est effectu^e pour tous Ies pixels de Timage 
NUMI, et on obtient Tlmage RFI dont la valeur de chaque pixel est egal au contraste 

20 C. 

L'avantage de ce filtre 1 est qu'il permet de selectionner des 
structures filiformes d'epaisseurs choisles et ainsi d'eliminer des structures plus 
epaisses. Par exemple, en imagerie medicale, ce filtre permet de selectionner.,une 
structure arborescente de vaisseaux sanguins et d'attenuer des structures plus 

25 epaisses comme une poche de sang. 

En outre, ce filtre peut effectuer une premiere selection permissive 
des pixels de I'image numerisee NUMI appartenant a une structure filiforme 
interessante, suivant un critere utilisant un parametre reglable X. Par exemple, Ies 
pixels selectionnes sont ceux dont la valeur est superieur a X fois la variance du 

30 bruit dans I'image fiitree par le dispositif decrit ci-dessus. Le parametre X peut etre 
choisi plus ou motns important selon le degre de selection desire. Le role de cette 
premiere selection est de limiter le nombre de pixels pris en compte dans les etf4pes 
ulterieures, afin d'eviter les problemes de surcharge de memoire de pile. 

Les figures 3a a 3d illustrent un exemple de technique destlnee a 

35 Testimatfon de la direction de chaque pixel d'image. La methode consiste a 

diagonaliser une matrice de Hessien HES, decrite sur la figure 3d, dont les elements 
diagonaux sont egaux a une premiere valeur Lxx et une deuxieme valeur Lyy et les 
autres elements sont egaux a une valeur Lxy. Pour chaque pixel de I'image 
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numerisee NUMI, on calcule une valeur Lxx, une valeur Lyy et une valeur Lxy. La 
figure 3a illustre le caicul de la valeur Lxx. Pour un pixel donne de Timage numerisee 
NUMI, on consldere un voisinage N55 de dimensions 5 pixels sur 5 pixels, centre sur 
ce pixel, et on effectue une convolution de ce voisinage N55 avec un masque Mxx. 
5 On obtient ainsi la valeur Lxx. La figure 3b illustre le caicul de la valeur Lyy. Pour le 
pixel donne de I'image numerisee NUMI, on effectue la convolution du voisinage 
N55 avec un masque Myy. On obtient ainsi la valeur Lyy. La figure 3c illustre le 
caicul de la valeur Lxy. Pour le pixel donne de I'image numerisee NUMI, on effectue 
la convolution du voisinage N55 avec un masque Mxy. On obtient ainsI la valeur Lxy. 

10 On effectue ensuite une diagonalisation de la matrlce de Hesslen HES et on retlent 
le vecteur propre U assocle a la plus petite valeur propre obtenue. Ce vecteur est le 
vecteur de direction tangentielle du pixel considere. On calcule ensuite la direction 
de ce vecteur, de telle sorte que Tangle obtenu est un angle compris entre 0 et 180 
degres. On divise ensuite cet angle par 10 et on arrondit le resultat obtenu a rentier 

15 superieur. On obtient ainsi un numero de secteur angulaire K pour le pixel 

considere. Dans cet exemple, il existe 18 numeros de secteur angulaire differents. II 
est possible d'utiliser un nombre different de numeros de secteur angulaire, par 
exemple en divisant Tangle du vecteur U par 5 au lieu de 10. 

La figure 4 illustre un exemple de dispositif d'analyse de connectivlte 

20 3 base sur la direction. L'exemple considere s'applique sur 6 pixels de Timage de 
direction DIRI numerotes de 1 a 6 (P1,P6), possedant 6 vecteurs de directions 
associes numerotes de 1 a 6 (U1,U6) et 6 numeros de secteur angulaire associes 
numerotes de 1 a 6 (K1,K6). Le principe de I'analyse de connectivlte est le suivant. 
Pour un pixel donne PN (N etant un nombre entler compris ici entre 1 et 6), on 

25 considere un voisinage NN s'etendant dans la direction du vecteur UN, sur une 
longueur I reglable. Dans le cas present, deux volsinages sont representes, 
symbolises par des ellipses en pointilles, Le voisinage Nl est un voisinage du pixel 
PI s'etendant dans la direction du vecteur Ul. Le voisinage N5 est un voisinage du 
pixel P5 s'etendant dans la direction du vecteur U5. Supposons que le pixel PI soit le 

30 premier pixel de Timage de direction DIRI traite par le dispositif d'analyse de 

connectivlte 3. La valeur de son numero de secteur angulaire Kl dans Timage de 
direction DIRI est comprise, par exemple, entre 0 et 18. Le pixel PI est alors 
labellise avec une etiquette egale a 19. Tous les pixels du voisinage Nl ayant un 
vecteur de direction associe proche de Ul seront labelltses avec cette meme 

35 etiquette 19, on dira quils appartiennent au label 19. Par exemple, on peut 

considerer qu'un vecteur UN d'un pixel PN du voisinage Nl est proche du vecteur Ul 
si les angles de ces 2 vecteurs different de moins de 10 degres, c'est a dire si le 
numero de secteur angulaire KN et le numero de secteur angulaire Kl different au 
maximum d'une unite. Dans Texemple considere, les pixels P2 et P3 appartiennent 
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au voisinage Nl. Le numero de secteur angulaire K2 etant egal au numero de 
secteur angulaire Kl et le numero de secteur angulaire K3 etant superleur de une 
unite au numero de secteur angulaire Kl, les pixels P2 et P3 appartiendront 
egalement au label 19. La meme operation est reallsee sur des voisinages des pixels 
5 P2, P3 et P4. II est alors facile de constater que les pixels P2, P3, P4, et P5 

appartiendront tous au label 19. Le pixel P6 possede un numero de secteur angulaire 
K6 superleur de 6 unites au numero de secteur angulaire K5. Par consequent, il 
appartiendra a un label different de 19, par exemple 20. L'operation est ainsi 
effectuee de fagon recursive sur tous les pixels de I'image de direction DIIU, jusqu'a 

10 ce que tous les pfxefs de I'image de direction DIRI solent labelllses. On notera 

qu'une fols qu'un pixel de Timage de direction DIRI a ete labellise, son Etiquette ne 
peut plus changer lorsque i'analyse de connectivite se poursuit. On obtient ainsi 
Timage CONI, constituee des differents labels auxquels appartiennent ies pixels. 

La derniere 6tape du procede consiste a selectionner des 

15 groupes de pixels dans I'image flltree RFI, afin de conserver les structures flliformes 
interessantes de cette image. Le princIpe de la selection de groupes de pixels est le 
suivant. On effectue les sommes \i des contrastes (dans I'image flltree RFI) de tous 
les pixels appartjenant au mSme label dans I'image connect^e CONI. On realise cette 
operation pour tous les labels de I'image CONI, et on considere le maximum Y de 

20 tous les \i calcules. On ne retlent que les pixels appartenant a un label tel que |i est 
superieur au quotient de Y par un parametre M reglabie. On obtient aInsI i'image 
traltee PRCI, constituee des pixels de I'image filtree RFI qui ont ete retenus api;es la 
selection de groupes de pixels. Le choix du parametre M permet d'obtenir differents 
degres de selection. En effet, si M est faible (en pratique de I'ordre de quelques 

25 unites), peu de groupes de pixels de I'image flltree RFI seront retenus dans I'image 
traitee PRCI. Si M est grand (par exemple de I'ordre de quelques centaines en 
pratique), la quasl-totalite des groupes de pixels de I'image filtree RFI sera retenue 
dans I'image traltee PRCL II faut noter que pour des structures flliformes epaisses 
comme une structure arborescente de yaisseaux sanguins, la selection de groupes 

30 de pixels donne de bons resultats du fait que les pixels voislns appartenant a une 
telle structure ont generalement la meme direction. 

La figure 5a illustre un systeme de traitement d'images destine a la 
mise en oeuvre de I'Invention. le systeme de traitement d'images comprend un 
disposltif 5 de traitement d'images selon I'Invention, une station de travail 6, un 

35 ecran 7, une souris 8 et un clavier 9. 

L'image numerisee NUMI est traitee par le dispositif 5 de traitement 
d'images, selon le procede decrit dans I'Invention. Ce dispositif 5 de traitement 
d'images peut etre implemente sous la forme d'un microprocesseur programmable 
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faisant partie de la station de travail 6 ou d'un circuit comprenant des elements 
comme des filtres, des memoires ou des operateurs logiques destines a realiser les 
differentes etapes decrites dans invention. L'ecran 7, la souris 8 et le clavier 9 
permettent la visualisation des images traitees et permettent le reglage des 
5 differents parametres regiables decrits dans Tinvention. Le systeme de traitement 
d'images peut etre relle h un dispositif de stockage afin de stocker les differentes 
Images obtenues au cours du traitement damages. 

La figure 5b lllustre un appareil d'imagerie medicale destine d la mise 
en oeuvre de I'lnvention. L'appareil comprend une source de rayons X 10, une table 

10 11 destin^e d recevoir un patient 12 a examiner, un systeme optlque 13 et un 
systeme de traitement d'images 14, 

Le systeme optlque 13 permet d'acquerlr des Images num^ris^s h 
partir des rayons X issus de la source de rayons X 10 ayant traverse le patient 12. Le 
systeme de traitement d'images 14 peut etre, par exemple, le systeme de traitement 

15 d'images decrit sur la figure 5a. 

La description ci-dessus en reference aux figures illustre Tinvention 
piutot qu'elle ne la limite. A cet egard, quelques remarques sont faites ci-dessous. 

II existe de nombreuses manieres pour selectionner des structures 
filiformes dans une image. La figure 2 illustre seulement une mise en oeuvre possible 

20 dans laquelle le filtre 1 utilise deux voisinages d'un pixel, avec des dimensions 
regiables. II est possible d'utiliser d'autres filtres pour selectionner des structures 
filiformes dans une image, prenant la forme de masques par exemple. 

II existe de nombreuses manieres pour effectuer une estimation de 
direction dans une image donnee. La figure 3 Illustre seulement une mise en oeuvre 

25 possible dans laquelle Testimatlon de la direction de chaque pixel s'appuie sur la 
diagonalisation d'une matrice de Hessien. II est possible d'utiliser d'autres 
techniques pour effectuer I'estimatton de la direction de chaque pixel, par exemple 
en utiiisant une vole discrete. 

II existe de nombreuses manidres pour effectuer une analyse de 

30 connectivity de pixels voislns bas^ sur leurs directions dans une image donnee. La 
figure 4 illustre seulement une mise en oeuvre possible dans laquelle I'analyse de 
connectivite de pixels voisins basee sur leurs directions utilise un voisinage d'un 
pixel s'etendant dans la direction du pixel, la taille de ce voisinage etant regiable 
grace a un parametre I. II est possible d'utiliser d'autres techniques pour effectuer 

35 I'analyse de connectivite basee sur la direction, avec des voisinages differents. 



Revendications 

1. Procede de traitement d'images dans des images contenant des 

5 structures filifortnes, comprenant en paralleie les etapes suivantes. 
une etape de filtrage desdites images ; 

une ^tape de decision destinee k selectionner les pixels d'image 
appartenant a une structure ffliforme interessante, ledit proced6 etant caracterise en 
ce que I'etape de decision comprend en paraJl^le une sous-etape d'estimation de la 

10 direction de chaque pixel d'image ainsi qu'une sous-etape d'analyse de connectivlte 
de pixels voisins basee sur leurs directions a la sortie de la sous-etape d'estimation 
de la direction de chaque pixel d'image et une sous-etape de selection de groupes 
de pixels, en fonction du resultat de la sous-etape d'analyse de conpectivite de 
pixels voisins basee sur leurs directions, a la sortie de I'etape de filtrage. 

15 2. Procede de traitement d'images selon la revendication 1, caract^ris^ 

en ce que ladite etape de filtrage desdites Images comprend une sous-etape de 
selection de pixels, les pixels d'une image seiectlonnes etant ceux dont le contraste 
est superieur a X fois la variance du bruit dans I'image, X etant un parametre 
reglabie par I'utilisateur. 

20 3- Procede de traitement d'images selon Tune des revendications 1 ou 2, 

caracterise en ce que ladite etape de filtrage desdites images utilise deux voislnages 
(Nl) et (N2) d'un pixel donne, I'ecart (GAP) entre ces voislnages etant reglabie par 
Tutilisateur. 

4. Procede de traitement d'images selon la revendication 3, caracterise 
25 en ce que la hauteur (H) et la longueur (L) desdits voislnages sont reglables par 

I'utilisateur. 

5. Procede de traitement d'images selon Tune des revendications 1 a 4, 
caracterise en ce que ladite sous-etape d'analyse de connectivity de pixels voisins 
basee sur leurs directions utilise un voisinage d'un pixel donne, ce volsinage 

30 s'etendant dans la direction du pixel considere, determinee lors de ladite sous-etape 
d'estimation de la direction de chaque pixel d'image. 

6. Procede de traitement d'images selon la revendication 5, caracterise 
en ce que la longueur dudit voisinage est reglabie par i'utilisateur. 

7. Procede de traitement d'images selon I'une des revendications 1 a 6, 
35 caracterise en ce que ladite sous-etape de selection de groupes de pixels utilise un 

parametre M reglabie par Tutilisateur, ce parametre M permettant de calculer la 
somme minimale des contrastes des pixels d'un groupe donne requise pour que ce 
groupe soit selectionne. 
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8. Programme d'ordinateur comprenant des instructions de code de 
programme pour ['execution des etapes du procede seion la revendicatlon 1 lorsque 
■edit programme est execute sur un ordlnateur . 

9. Systeme de traitement d'images comprenant un ordinateur destine a 
la mise en oeuvre d'un programme d'ordinateur seion la revendicatlon 8 ou un circuit 
destine a ia mise en ceuvre du procede de traitement d'images seion I'une des 
revendications 1 a 7, un dispositif de projection des images traltees seion ledit 
procede et eventuellement un dispositif de stockage desdites images. 

10. Appareil d'imagerie medlcale comprenant un dispositif d'acquisition 
damages et un systeme de traitement d'images seion la revendicatlon 9. 



3/5 



:3 





